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The scattering theory approach makes it possible to carry out exact calculations of Casimir en- 
ergies in any geometry for which the scattering T-matrix and a partial wave expansion of the free 
Green's function are available. We implement this program for the case of a perfectly conducting 
elliptic cylinder, thereby completing the set of geometries where electromagnetic scattering is sepa- 
rable. Particular emphasis is placed on the case of zero radius, where the elliptic cylinder reduces 
to a strip. 
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I. INTRODUCTION 

Formulating the Casimir energy in terms of scattering theory has made it possible to efficiently reduce quantum field 
theory calculations to standard problems in quantum mechanics and electromagnetism. By expressing the "TGTG" 
form of the Casimir energy [l[ in appropriate scattering bases, one can calculate the Casimir interaction energy of 
a collection of objects as a combination of the objects' scattering amplitudes (T-matrices) together with universal 
translation matrices, which are obtained from a mode expansion of the free Green's function [2|-y. The former are 
computed for each object individually, while the latter depend only on the objects' relative positions and orientations. 
As a result, the Casimir energy can be computed for any collection of objects for which the scattering T-matrix is 
available within a standard scattering basis. This approach allows for exact calculations, extending earlier results 
using asymptotic expansions 0, @ and results from scalar theories 7 9]. It can also be applied in the weak coupling 
approximation [To| . For objects without special symmetries, however, one must ultimately turn to computational 
methods to compute either the T-matrix or the associated Green's functions flll - flij . 

With sufficient symmetry, the exact T-matrix can take an analytically calculable form, greatly reducing the amount 
of computation required. This reduction has made it possible to apply the scattering method to efficient computations 
of the Casimir energy for planes [l5| , spheres and ordinary cylinders @, 0, [l|| [TtJ , parabolic cylinders [H, EH , and 
wedges and cones [2fJ. Here we complete the set of separable geometries in electromagnetism by treating the case 
of an elliptic cylinder. This geometry has been investigated for microfabricated materials using a Lifshitz formula 
approach in Ref. [2l[ and has been used to study Casimir self-energies in Refs. [22l [23j. 

II. SCATTERING IN ELLIPTIC CYLINDER COORDINATES 

We begin by formulating scattering theory in elliptic cylinder coordinates, 

x = d cosh [i cos 9 

y = d sinh fi sin 9 , (1) 

where 2d is the interfocal separation of our elliptic cylinder coordinates, 9 is the analog of the angle in ordinary 
cylindrical coordinates, and /i is the analog of the radius, with 



r=Vx 2 +y 2 = d\l ^ >-e" (2) 



' cosh 2(i + cos 29 d „ 
> — v *" 1 .y — ■* \i ~ 

as /i — > oo. 

We use separation of variables to form solutions of the Hclmholtz equation — V 2 V>( r ) = k 2 ip(r) as products of 
functions of (i, 9 and z individually. For the functions of z, we have ordinary complex exponentials e lk ' z , which will 
multiply angular functions of 9 and radial functions of \x. Since we have parity symmetry, we can choose our angular 
solutions to be either even or odd under reflection across the x-axis, 9 — > —9. Unlike the ordinary cylinder case, the 
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elliptic angular solutions depend on the wave number k, and the elliptic radial solutions associated with the the odd 
and even angular solutions differ and depend on the wave number and radius separately, rather than only on the 
product kr. For q = ^-(k 2 — k 2 ), the angular solutions are the even and odd angular Mathieu functions ce m (9, q) and 
se m (9,q), which are the analogs of cosmO and sinmf? respectively. As in the case of ordinary cylindrical coordinates, 
for the even functions m runs from to oo, while for the odd functions m runs from 1 to oo. For the corresponding 
radial functions, we have both the even and odd first kind solutions Je m (n, q) and Jo m (n, q), the analogs of the Bcssel 
function J m {\/k 2 — k 2 r), and the even and odd outgoing wave solutions He m (n, q) and Ho m (n, q), the analogs of the 

Hankel function Hm\\/k 2 — k 2 r). We will normalize the Mathieu functions so that they obey the same orthonormality 
conditions as their cylindrical analogs, except that the m = even angular function will be normalized so that its 
root mean square average value is l/v2 (the same as for all the other angular functions) instead of cosO = 1. As a 
result, we have 



2jv 



ce m (9,q) 2 d9 



2jv 



se m (0, qfdO 



(3) 



with the radial functions normalized to coincide with their cylindrical analogs asymptotically. Our notation and 
normalization match that of Ref. [24| , which defines Mathieu functions following the conventions of Abramowitz and 
Stegun [25l |. but uses a modified notation that is more closely analogous to the ordinary cylinder case. We will make 
use of identities for elliptic cylinder functions found in standard references [25l - f27l ] . 
The key ingredients for our calculation will be the free Green's function 



G( ri ,r 2 ,k) = / -f- 



° dk z i 
2^2 



oo 



m=0 



ce m (0i,q)ce m (9 2 , q)Je m (n<,q)He m (n>,q) 



E 

m—1 



- se m (9 l7 q)se m (9 2 ,q)Jo m (n<,q)Ho m (n >7 q) 



where /i< (n>) is the smaller (larger) of ni and fx 2 , and the expansion of a plane wave, 

oo oo 

2 ^2 i m ce m (4>, q)ce m (9, q)Je m (n, q) + 2 ^ i m se m ((f),q)se m (9,q)Jo m (p,q) 



Ah-r 



m— 



m—1 



(4) 



(5) 



where [A, 9, and z are the elliptic cylinder coordinates of r and <j) : 
xy-plane, with k 2 = k 2 + k 2 + k 2 . 

We will work on the imaginary A:-axis k = in, so that k y = iyj k 2 



arctan -p- is the angle of k = {k x , k y ,k z ) in the 



k' 2 + k 2 and q = —d 2 (K 2 + k 2 )/4 is negative. As 



a result, it is convenient to rewrite these expressions in terms of modified radial functions, 

p0C dk z 1 

27T 7T 



G(n,r 2 , k) 



OO 
OO 

■E 

m—1 



^ ce m (6 1 ,q)ce m (9 2 ,q)Ie m (n < ,-q)Ke m (fJ, > ,-q) 



xn— 



,q)se m (9 2 ,q)Io m (n < ,-q)Ko m (n > , -q) 



(6) 



and 



Ak-r 



Ah z z 



oo 

7Ti=0 



{-l) m ce m ((t>, q)ce m (9, q)Ie m (n, -q) + 2 ^ (-l) m se m ((j), q)se m (9, q)Io m (fi, -q) 



(7) 



where Ie m (/j,,-q) = i m Je m (n,q), Io m {^,-q) = i m Jo m (^,q), Ke m (p,-q) = i rn+1 ^He m (n,q) and Ko m (/j,,-q) 



■m+l 7r 



Ho m (/i, q) are the modified outgoing radial functions. 



We will consider scattering with Dirichlct and Neumann boundary conditions on an elliptic cylinder of radius /io. 

■e,o 

mk z m' k' 



For the scattering amplitudes, we have , k , = 2ir5(k z — k' z )S mm 'T^° , with 



' 771 

T 

' ID 



Ie m (no, -q) 
Ke m (no,-q) 

te'm (mo, -q) 



T" 
T" 



Ke' rn (no,-q) 

where prime indicates a derivative with respect to ji 



Iom (no, -q) 

Ko m (no,~q) 

/o m(M0:-g) 
K °'rn (V0,-q) 



(Dirichlet) 
(Neumann) 



(8) 
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III. ELLIPTIC CYLINDER AND PLANE 



To consider the elliptic cylinder's interaction with a plane, we will need to connect the elliptic cylinder and planar 
geometries. To do so, we make use of the expression for the free Green's function in Cartesian coordinates for y- 2 > y%, 



G(r 1 ,r 2l k) = 



dk z 
2^' 



,ife z (z 2 — zi) . 



4-7T 



dk 



HL e i { k x( x 2-xi)+k y (y 2 -yi)) 



(9) 



where k y = y 'k 2 — k 2 — k 2 = i^J ' k 2 + k 2 + k 2 . We equate Eq. © to the Green's function in Eq. ([5]), expand the 
plane wave e lk r2 in Eq. ([9]) using Eq. 0, make the substitution k x — > —k x , and finally use the orthogonality of the 
regular elliptic cylinder solutions to equate both sides term by term in the sums over m. The result is an expansion 
for the elliptic outgoing wave solutions in terms of plane waves for y < [25| , 



ce m {6 1 q)Ke m { r i 1 -q)e lk ~~ z 
se m (9,q)Ko m {n,-q)e tk * z 



dk x 



2k 



■ce m (<t>, q) 



g—ik y y-\-ik x x ^Lk z z 



(10) 



The quantities in brackets represent the translation matrix elements, which we must then multiply by the normalization 

£Tcllipti< 



factor 



where we can read off C" 



= A /i and CV pla 



= . from the expressions for the free Green's 



function in Eqs. ((6]) and (|9|). 

Finally, the T-matrix elements for the plane in Cartesian coordinates are simply T p = ±1 for Neumann and 
Dirichlet boundary conditions respectively. (For more general boundary conditions on the plane, this scattering 
amplitude would be a function of k x .) 

We have now obtained the T-matrix elements, which describe how waves scatter off each object individually, and 
the translation matrix elements, which convert the scattering bases between the two objects. As a result, we are 
prepared to assemble these ingredients into the result for the full Casimir interaction energy per unit length. We 
consider a perfectly conducting plane oriented perpendicular to the y-axis and a perfectly conducting elliptic cylinder 
with its z-axis parallel to the plane, its center at a distance H from the plane, and its major axis at an angle ip to 
the plane, as shown in Fig. [TJ This angle represents a rotation of the elliptic cylinder coordinates and fi relative to 
the Cartesian coordinates x and y, which we then implement in Eq. (|10[) by adding a constant shift ip to the angle 
d> = arctan in the translation matrix elements. 




H 

V 



FIG. 1: Geometry for the elliptic cylinder and plane. 

For a particular choice of boundary conditions, we can now use the approach of Rcfs. [2- 4] to write the Casimir 
energy per unit length as 

e_ = p g p logdct ( ^ _ ^ i ipy^TPO* ) , 



hcL J n 2n J m 2ir 



where the matrix determinant runs over x, x' = o,e with m = 0, 1, 2, 3 . . . for \ = e and m = 1,2,3... for x = o, and 
similarly for m! and x' ■ The translation matrices and reverse translation matrices are given by 

K*. = ce m (cb + ip, q) e^ H W mkx = ce m (-0 + <p, q) e? h * H (12) 
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for the even modes and 

U° mkx = se m {4> + tp, q) e^ H U° mhx = se m (-0 + <p, q) ^ H (13) 

for the odd modes. 

We can then change the integration variable from k x to u = i (0 — ^) and combine the k and k z equations into a 

I rl 2 2 

single integral over p = k 2 + fc|, so that q = 2-. We obtain 



4z = i 1 



(i 



vsL - vkr p r due-**«*» c y (l+ iu+ *> q ) IT 



(14) 



where we choose ce m for x = e an d se m for x = °i an d similarly for m! and x'. The full electromagnetic Casimir 
energy is the sum of this result for Dirichlet conditions on both surfaces and for Neumann conditions on both 
surfaces. Note that the established result for an ordinary cylinder [l6[ can be obtained from this expression by 
replacing the elliptic functions with their ordinary cylindrical analogs, combining the even and odd modes using 
cosh mu cosh m! u + sinh mu sinh m'u = cosh(m + m')u, and employing the integral identity 

poo 

K n {a)= e- acoshu coshnudu. (15) 



There are several special cases of interest in which the calculation simplifies: 

• Plane perpendicular to the ellipse's major axis. 

For (p = 7r/2, the elliptic cylinder's major axis runs perpendicular to the plane. By the reflection symmetry 
across the y-axis, the even and odd sectors decouple, and we can compute the Casimir energy by considering 
the odd and even elliptic modes separately. 

• Plane parallel to the ellipse's major axis. 

For ip = 0, the elliptic cylinder's major axis lies parallel to the plane. This case also has reflection symmetry 
across the y-axis, but this symmetry does not correspond directly to the symmetry of the even and odd Mathieu 
functions. Instead, the even Mathieu functions of even order and the odd Mathieu functions of odd order are 
symmetric under this transformation, while the odd Mathieu functions of even order and the even Mathieu 
functions of odd order are antisymmetric. (This is the same symmetry structure as the ordinary trigonometric 
functions have when their argument is displaced by 7r/2.) Thus we can again decompose the problem into two 
independent sectors, consisting of the modes for which the parity of the elliptic functions matches the parity of 
m, and the modes for which they are opposite. 

• Zero radius cylinder. 

An elliptic cylinder with fiQ = becomes a strip of width 2d, allowing us to study the effects of edges (]~8l — 
I20I l28l . |29| . In that case we have 7J£ = for a Dirichlet boundary and 7^ = for a Neumann boundary, since 
in these cases the free modes already obey the boundary condition at the surface. These modes therefore give 
zero contribution to the Casimir energy in this case, and can be omitted from the calculation. 



IV. NUMERICAL RESULTS 



We can now compute the Casimir energy by straightforward numerical integration of Eq. (1141). To compute the 
modified radial functions needed for the scattering amplitude, we use the package of Alhargan [30L l3l| . (Standard 
packages such as Maple and Mathcmatica only implement the angular Mathieu functions of the first kind. Although 
the radial functions are related to the angular functions with imaginary argument, without an implementation of the 
second kind angular function we cannot take advantage of this relationship to compute the functions needed to for 
the scattering amplitude.) The angular functions arising from the translation matrix, on the other hand, need to be 
computed for complex arguments, which are not supported directly in the Alhargan package. Fortunately, since only 
the first kind angular functions are required, we can use the implementation in Mathcmatica, which supports fully 
complex arguments. As a final complication, because of problems with the Mathieu function routines in the current 
version of Mathcmatica for the case where the parameter q < 0, we make use of the identities 

f (— l)^ce m (j - fx, -q) for m even 
ce m {fj.,q) = < 'nd v % \ t jj 

L (— 1) 2 se m [t£ — //, — q) for m odd 
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FIG. 2: Electromagnetic Casimir interaction energy for a perfectly conducting strip opposite a perfectly conducting plane, as 
a function of the orientation angle <p. The distance H from the center of the strip to the plane is twice the distance from the 
center of the strip to the edge of the strip, H — 2d. The solid line shows the proximity force approximation. 



f (-1)' 2 " 1 se m (f - n,-q) for m even 
[ (-1) 2 ce m (-| — /i, —q) lor m odd 

so that we only need to compute the angular functions for — q > 0. As a result, after importing the Alhargan routines 
for the modified radial functions, it is possible to carry out the full calculation within Mathcmatica. Because of 
limitations in the ability of the angular routines to handle large imaginary arguments, however, it was not possible 
to extend the calculation to very small separations. 

Figure [2] shows the orientation dependence of the Casimir interaction energy for a perfectly conducting strip (an 
elliptic cylinder of zero radius) for the case where the distance H from the center of the strip to the plane is twice 
the distance from the center of the strip to the edge of the strip, H = 2d. Because higher values did not change the 
results appreciably, the matrix determinants were truncated at m max = 8. We see that the lowest energy occurs for 
if = 7r/2, when the strip is perpendicular to the plane. As expected, the result for the energy per unit length in this 
case, i^r = —0.00637, is less negative than the —0.00674 one finds f° r the case where the strip is extended to 

an infinite half-plane whose edge maintains the same distance H — d = d from the infinite plane. We note, however, 
that if we subtract the contribution from a half-plane at distance H + d = 3d from the result for the half-plane at 
distance H — d = d to account for the missing remainder of the half plane, we obtain —0.00674 ■ | = —0.00599, which 
underestimates the magnitude of the true result for the strip. We also compare these results to the proximity force 
approximation (PFA), 



4°]m _JP_ r dcosv dx 7T 2 Hd cosy) 

HcL " 720 J_ d cos v (H + xtanf) 3 ~ 360 (#2 _ <p sin 2 ^ 



(17) 



which gives a good approximation for tp = but goes to zero at tp = ir/2. For ip ^ the derivative expansion 
correction to the PFA [32|, [33| is also invalid, because of the sharp curvature at the point of closest approach. 

Figure [3] shows the Casimir interaction energy for a strip oriented parallel to a plane as a function of the distance 
to the plane. The energy is shown as a ratio with the PFA result (in this case the correction from the derivative 
expansion vanishes). As in the case of the ordinary cylinder [l6|, the PFA is an underestimate at large distances, but 
at short distances the exact result approaches the PFA result from below. These calculations were carried out with 
the matrices truncated at several different values of m max up to m max = 16, with the final result then obtained by 
extrapolating these results for m max — > oo. 



V. DISCUSSION 



We have computed the Casimir interaction energy for an elliptic cylinder, the last remaining geometry for which 
electromagnetic scattering is separable. For a plane, cylinder, and sphere, the problem remains separable even for 
a dielectric, while for a parabolic cylinder, elliptic cylinder, wedge, and cone only perfect conductors can be solved 
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FIG. 3: Ratio of the electromagnetic Casimir interaction energy to the proximity force approximation (PFA) for a perfectly 
conducting strip of width 2d parallel to a perfectly conducting plane, as a function of the separation H. As in the case of an 
ordinary cylinder [l(|, the ratio is a nonmonotonic function of H. 



exactly. However, the scattering method is particularly useful in these latter cases, because they contain sharp limits 
in which the PFA is invalid. In principle, it should be possible to extend the elliptic cylinder result to a hyperbolic 
cylinder in the same way as the wedge is obtained from the ordinary cylinder and the cone is obtained from the sphere, 
but at present there do not appear to be routines available for computing all the Mathieu functions of complex order 
that would be needed for such a calculation. 

Focusing on the limit in which the elliptic cylinder becomes a strip has made it possible to study the orientation de- 
pendence of the Casimir force, to show how the PFA depends on distance and angle, and to observe non-superposition 
effects in the perpendicular configuration (where the PFA is invalid). With improvements to the available routines 
for computing Mathieu functions, this calculation could offer an independent check of the edge correction that was 
obtained for half-planes [HI, Ell EH • More generally, this calculation establishes another addition to the toolbox of 
Casimir problems that can be cast into analytically tractable form. 
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